Modulational instability criteria for two-component Bose— Einstein condensates * 
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The stability of colliding Bose-Einstein condensates is investigated. A set of coupled Gross- 
Pitaevskii equations is thus considered, and analyzed via a perturbative approach. No assumption 
is made on the signs (or magnitudes) of the relevant parameters like the scattering lengths and 
the coupling coefficients. The formalism is therefore valid for asymmetric as well as symmetric 
coupled condensate wave states. A new set of explicit criteria is derived and analyzed. An extended 
instability region, in addition to an enhanced instability growth rate is predicted for unstable two 
component bosons, as compared to the individual (uncoupled) state. 
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I. INTRODUCTION 

Bose-Einstein condensation of dilute gases in traps has 
attracted a great deal of interest recently, as witnessed 
in recent reviews and monographs Q, ■ Mean-field the- 
ory provides a consistent framework for the modeling of 
the principal characteristics of condensation and eluci- 
dates the role of the interactions between the particles. 
A generic theoretical model widely employed involves the 
Gross-Pitaevskii equation, which bears the form of a non- 
linear Schrodinger-type equation, taking into account bo- 
son interactions (related to a scattering length a) , in ad- 
dition to the confinement potential imposed on the Bose- 
Einstein condensates (BECs) in a potential trap. The 
scattering length a, although initially taken to be posi- 
tive (accounting for repulsive interactions and prescrib- 
ing condensate stability) , has later been sign- inverted to 
negative (attractive interaction) via Feshbach resonance, 
in appropriately designed experiments. This allowed for 
the prediction of BEC state instability, eventually leading 
to wave collapse, which is only possible in the attractive 
case (a < 0) [lj . As expected from previous know-how on 
problems modelled by generic nonlinear Schrodinger-type 
equations (in one or more dimensions), the analysis of 
BEC dynamics revealed the possibility for the existence 
of collective excitations including bright- (for a < 0) 
and dark- (holes, for a > 0) type envelope excitations, 
as well as vortices, which were quite recently observed 
in laboratories 0, 0, IE 0- The evolution of coupled 
("colliding") BEC wavepackets was recently considered 
in theoretical and experimental investigations 0, 0> @ • 
Pairs of nonlinearly coupled BECs are thus modeled via 
coupled Gross-Pitaevskii equations, involving extra cou- 



* Submitted to European Physics Journal B. 
^Electronic address: ioannis0tp4.rub.de ; 
http : //www . tp4 . rub . de/~ioannis. 

* Electronic address: ps@tp4.rub.de ; 
http : / /www . tp4 . rub . de/ ^ps . 



pling terms whose sign and/or magnitude are a priori not 
prescribed. Although theoretical modeling, quite natu- 
rally, first involved symmetric pairs of (identical) BECs, 
for simplicity, evidence from experiments suggests that 
asymmetric boson pairs deserve attention 

In this paper, we investigate the stability of a nonlin- 
early coupled BEC pair, from first principles. Both BECs 
are assumed to lie in the ground state, for simplicity, al- 
though no other assumption is made on the sign and/or 
magnitude of relevant physical parameters. We shall de- 
rive a set of general criteria for the stability of BEC pairs 
(allowing for asymmetry in the wave functions). 



II. THE FORMALISM 

The wave- functions ipi and ^2 of two nonlinearly in- 
teracting BECs evolve according to the coupled Gross- 
Pitaevskii equations (CGPEs) 



dt ' 2rni 



+/xiV>i = , 



dt 2m,2 



(1) 



where V 2 = d 2 /dx 2 + d 2 /dy 2 + d 2 /dz 2 is the Laplace 
operator (a three-dimensional Cartesian geometry is con- 
sidered, for clarity). Here rrij represents the mass of 
the jth condensate. According to standard theory, the 
nonlinearity coefficients Vjj are proportional to the scat- 
tering lengths Oj via Vjj = Airhaj /rrij, while the cou- 
pling coefficients Vji are related to the mutual interac- 
tion scattering lengths aji via Vji — 2Trhaji/mji, where 
rriji = rrijmi/ (rrij +mi) is the reduced mass. The (linear) 
last terms in each equation involve the chemical potential 
fij, which corresponds to a ground state of the conden- 
sate, in a simplified model. These terms may readily be 
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eliminated via a simple phase-shift transformation, viz. 
ipj = ipj exp(ifijt) (j = 1, 2); this is however deliberately 
not done at this stage, for generality. Nevertheless, one 
therefore intuitively expects no major influence of the 
chemical potentials on the coupled BEC dynamics (at 
least for the physical problem studied here). 

III. LINEAR STABILITY ANALYSIS 



where the coupling is expressed via fl 4 = Q 2 2 Q 21 = 
M 12 M 2 i/h 4: in the right-hand side of Eq. (4). We stress 
that this dispersion relation (which is independent of the 
chemical potentials fij) relies on absolutely no assump- 
tion on the sign or the magnitude of rrij, Vjj and Vji. 

IV. MODULATION AL INSTABILITY OF 
INDIVIDUAL BECS 



We shall seek an equilibrium state in the form ipj — 
ipjoexp[iipj(t)], where tpjo is a (constant real) amplitude 
and <pj(t) is a (real) phase, into the CGP Eqs. @J. 
We then find a monochromatic (fixed- frequency) Stokes' 
wave solution in the form: ifj(t) = Qjot, where 



jo 



fij 



iorj ^1 = 1,2. 

Let us consider a small perturbation around the sta- 
tionary state defined above by taking ipj — (ipjQ + 
eipji) exp[i(pj(t)}, where ipji(v, t) is a complex number de- 
noting the small (e <C 1) perturbation of the slowly vary- 
ing modulated bosonic wave-functions (it includes both 
amplitude and phase corrections), and ipj(t) is the phasor 
defined above. Substituting into Eqs. 0) and separating 
into real and imaginary parts by writing tpji = a,j + ibj, 
the first order terms (in e) yield 



In the vanishing coupling limit, i.e. for Vji — > 0, the 
dispersion relation J3J gives fi&,± = ±fij (j — 1,2). 
Absolute stability is ensured if Vjj > 0. On the other 
hand, if Vjj < 0, a purely growing unstable mode oc- 
curs (viz. f2 2 < 0) for wavenumbers below a critical 
value kj tCr — 2(mj\Vjj\) 1 ' 2 \tpjo\/ft. The growth rate a — 
iy— f2? attains a maximum value a max = \Vjj\\ipjo\ /ft 
dbti k - — - kj ct / 

Recalling the definitions of Vjj, we see that a repul- 
sive/attractive scattering length (i.e. positive/negative 
Vjj) prescribes a stable/unstable (single) BEC behavior. 
In the following, we shall see how this simple criterion 
for stability (Vjj > 0) is modified by the presence of in- 
teraction between two he condensates. 



V. MODULATIONAL INSTABILITY OF 
COUPLED BECS 



dbj h 2 2 2 
- ft-^j- + a i ~ 2V ]j^w a j - ZVji-ipjoipmai = , 

where j and l(^= j) = 1,2 (this will be henceforth un- 
derstood unless otherwise stated). Eliminating bj, these 
equations yield 



2 3 2 



+ 



ft 2 ( h 2 



dt 2 2m\ \2mi 



Oi 



The dispersion relation Q takes the form of a bi- 
quadratic polynomial equation 



n 4 k - ra 2 + d = o, 

s 2 _ 



(5) 



4 _ 



where T = TrM/?i = Of + O 2 . and D = DctM/fr 

are related to the trace and the deter- 
minant, respectively, of the matrix M. Eq. J5j has the 
solution 



i 



[T±(T 2 -4£>) 1/2 ] 



(6) 



or 



V 12 \iP w \\ip 20 \V 2 a 2 =0, 

mi 



(3) 



(together with a symmetric equation, obtained by per- 
muting 1 «-> 2). We now let aj — aj exp[i(k • r — 57^)] + 
complex conjugate, where k and f2fc are the wavevec- 
tor and the frequency of the modulation, respectively, 
viz. d I dt — ► — iVLk and d/dx n — > ik n (x n = {x,y,z} 
for n = 1,2,3) i.e. d 2 /dt 2 -> -Q 2 , and V 2 -> -fc 2 . 
After some algebra, we obtain the eigenvalue problem: 
Ma = (hw) 2 a., where a = (ai,a 2 ) T , and the matrix ele- 
ments are given by Mjj = ej(ej+2Vjj\ipjo\ 2 ) = ft 2 ^ 2 and 
Mji = — 2ejVji\ipjo\ \ipio | = ft 2 Q 2 i\ where we have defined 

ej = Ti 2 k 2 /2mj. The frequency uo and the wave number 
k are therefore related by the dispersion relation 



n 2 ) (n 2 - n 2 ) = d* 



nl ± = \(n 2 + n 2 ) ± \ [(n 2 - n 2 ) 2 + 4^] 1/2 . (7) 

We note that the right-hand side is real/complex if the 
discriminant quantity A = T 2 — AD is positive/negative, 
respectively. 

Stability is ensured (for any wavenumber k) if (and 
only if) both solutions fl 2 , ± are positive. This is tanta- 
mount to the following requirements being satisfied si- 
multaneously: T > 0, D > and A > 0. Since the three 
quantities T, D and A are all even order polynomials 
of k, one has to investigate three distinct polynomial in- 
equalities. The stepstones of the analysis will be outlined 
in the following, though trying to avoid burdening the 
presentation with unnecessary details. 

First, the sign of T = k 2 [(h 2 k 2 /4)^.(l/m 2 ) + 
(4) J2j Vjjl^jol 2 / m j] ( see definitions above) depends on (the 
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sign of) the quantity J2jVjj\i ; jo\ 2 /' m j which has to be 
positive for all k, in order for stability to be ensured (for 
any ipjo and k). This requires that 

Vu > and V 22 > . (8) 

Otherwise, T becomes negative (viz. ft? _ < 0, at 
least) for k below a critical value k cr ,i = V Kx, where 
Kx = 4(-E J -^il^o| 2 M)/[ft 2 Ei(V^)] > (cf. the 
single BEC criterion above); this is always possible for 
a sufficiently large perturbation amplitude IV'iol if, say, 
Vrx < (even if V22 > 0). Therefore, only a pair of two 
repulsive type BECs can be stable; the presence of one 
attractive BEC may de-stabilize its counterpart (even if 
the latter would be individually stable). 

Second, D — Clf Of— fif^li is an 8th-order polynomial 
in k, which can be factorized as D ~ fc 4 (fc 4 + bk 2 + c), 
where b — ^^,j{nijVjj)/h and c = 16mim 2 (VnV 22 — 

^12^21) |^io| 2 1^20 17^ (note that b 2 - 4c > 0). The 
stability requirements b > and c > (in order for 
D to be positive for any value of k > 0) amount to 
miVn + JTI2V22 > and 

^11^22-^12^21 >0, (9) 

respectively. Only the latter condition for stability has to 
be retained, since the former one is automatically covered 
by (8) above. To be specific, solving D = for k 2 = K 2i ±, 
viz. K 2 ,± = [-b ± {b 2 - 4c) 1 / 2 ]/2, we see that: 

(i) if b < < c, then < K 2 - < K 2 ,+, and D < for 
y/K2~ < k < y/K2,+ (instability for short wavelengths); 

(ii) if c < (regardless of 6), then -K2,- < < K 2y+ , and 
D < for < k < y/K^; 

(iii) if b > and c > 0, then _ff 2 ,- < ^2,+ < 0, so that 
D > 0. 

We see that this kind of instability, - i.e. if the criterion 
(9) is not met, is due to the mutual interaction potential 
Vji among the bosons. 

Finally, the positivity of A = T 2 - AD = (ft 2 - VL 2 ) 2 + 
4^2^21 is only ensured (for every value of k and |t/>jo|) 
if ft? 2 fi 2i ~ ^12^21 > 0, i.e. if 

V 12 V 2 i>0. (10) 

If this condition is not met, the solution © above has a 
finite imaginary part, which accounts for amplitude in- 
stability due to the external perturbation. For rigour, 
we note that A bears the form A = fc 4 (c4fc 4 — c 2 k 2 + cq) 
(where C4 > 0; the complex expressions for c n are omit- 
ted). If A' ee c\ - 4c c 4 F12F21 < 0, i.e. if (10) is 

met, then A > for any value of k. If A' > 0, on the 
other hand, denoting K 3t ± = [c 2 ± (c 2 — 4C0C4) 1 / 2 ]/ (2gi), 
we find that: 

(i) stability is only ensured (since ^3,- < ^3,4- < < k 2 ) 
if c 2 < < c (nevertheless, this condition depends on 



the perturbation amplitudes \ipjo\ and may always be vi- 
olated) . 

(ii) Again, a finite unstable wavenumber interval k S 
{y/K$-, y/K$, + ) is obtained for C2 > and Co > 0. 

(iii) Finally, instability will be observed for k £ 
(0, ^JK 3 ,+) if c < (regardless of c 2 ). 



VI. CONCLUSIONS 

Summarizing, we have derived a set of explicit criteria, 
(8) to (10) above, which should all be satisfied in order 
for a boson pair to be stable. Therefore, an interacting 
BEC pair is stable only if the interaction potentials sat- 
isfy V11 > and V 22 > and VnV 22 > V 12 V 2 i > 0. If 
one criterion is not met, then the perturbation frequency 
develops a finite imaginary part and the solution blows 
up in time. A few comments and qualitative conclusions 
should however be mentioned. 

First, for a symmetric stable boson pair, viz. Vu = 
V22 > and V12 = V21, stability is ensured if V 2 2 < V 2 X . 
Second, if one BEC satisfies Vjj < 0, the pair will be 
unstable: only pairs consisting of stable bosons can be 
stable. Interestingly, in the case of individually unstable 
BECs (viz. Vjj < 0, for j =1 or 2), the instability charac- 
teristics are strongly modified. For instance, in the case 
of a symmetric unstable boson pair (viz. Vu = V 22 < 
and mi = m 2 ), an extended unstable wavenumber region 
and an enhanced growth rate can be obtained, as can be 
checked via a tedious calculation; cf. Fig. 2] Further- 
more, we have pointed out the appearance of secondary 
instability "windows" , i.e. unstable wave number inter- 
vals beyond (k cr , k' cr ), where k cr 7^ 0. 

These results above follow from a set of explicit stabil- 
ity criteria. Both BECs were assumed to lie in the ground 
state, for simplicity, although no other assumption was 
made on the sign and/or magnitude of the relevant phys- 
ical parameters. Naturally, a future extension of this 
work should consider the external confinement potential, 
imposed on the trapped condensates. Our results can 
be tested, and can hopefully be confirmed, by designed 
experiments. 
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FIG. 1: The growth rate 7 versus the wavenumber k (in units 
of I Vjj I \tpjfi\ 2 /h and 2irij\Vjj \ \ipj,o\/fi, respectively) for a 
symmetric pair of coupled unstable (Vji < for j,l — 1,2) 
BECs (upper curve) as compared to the single BEC case 
(lower curve). 



